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EIGENVALUES OF SEBA BILLIARDS WITH LOCALIZATION 
OF LOW-ENERGY EIGENFUNCTIONS 

MINJAE LEE 


Abstract. We study the localization of eigenfunctions produced by a point 
scatterer on a thin rectangle. We find an explicit set of eigenfunctions localized 
to part of the rectangle by showing that the one-dimensional Schrodinger oper¬ 
ator with a Dirac delta potential asymptotically governs the spectral properties 
of the two-dimensional point scatterer. We also find the rate of localization in 
terms of the aspect ratio of the rectangle. In addition, we present numerical 
results regarding the asymptotic behavior of the localization. 


1. Introduction 

The Seba billiard, a point scatterer on a two-dimensional plate, was first intro¬ 
duced in ^ as a simple model to study quantum chaos. It is formally given by the 
Schrodinger operator 

(1.1) -A-|-c(5xo 

where A is the Dirichlet Laplacian, Xq is a point in the plate and c S M but requires 
renormalization or spectral theory to be properly defined. More precisely, a point 
scatterer at Xq is defined as a self-adjoint extension [5] of the Laplacian whose 
domain consists of functions vanishing near Xq. Such an extension is not unique 
and we can parametrize all possible self-adjoint extensions with a single parameter 
a € (— 00 , 00 ] which can be considered as a number related to the strength of the 
point scatterer. See miisiis] for further developments. 

In this paper we consider a point scatterer on a thin rectangle. The spectral 
theory of the Laplacian on a thin domain has been studied in mathematical physics 
as a model for waveguides, nanotubes, quantum wires and integrated circuits. Such 
objects are called quasi one-dimensional since its length in one direction is signif¬ 
icantly larger than those in the other directions. The simplest example of such 
domain is a rectangle in which the eccentricity, the ratio of the width to the height, 
is very large. There are other examples such as a cylindrical shell, e-neighborhood 
of a curve and e-neighborhood of a graph. Then it is natural to ask if a ID model 
on a curve or a graph corresponding to the thin domain can approximate the spec¬ 
tral properties of the 2D operator such as the Laplacian. This question for the 
Laplacian has been answered in various contexts. See [7] for details. See also 
[n m [m nn nzi for the spectral properties of the Laplacian on other kinds of 
domains. 

In [5] , we showed numerically that a point scatterer located inside a thin rectangle 
acts as a barrier so that the low-energy eigenfunctions except the lowest one get 
localized as the eccentricity tends to infinity. This work was motivated by the 
numerical study of Filoche and Mayboroda [3] for bi-Laplacian A^ on a rectangle 
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with a point removed. This fourth order operator exhibits localized modes under the 
boundary conditions requiring the mode and its gradient to vanish at the boundary. 

In this paper, we provide an explicit set of eigenvalues and corresponding strength 
parameters of the point scatterer exhibiting the localization of eigenfunctions with 
the sharp estimate of the error term. Interestingly, we observed that for low-energy 
eigenvalues and eigenfunctions there exists asymptotic correspondence between the 
Seba billiard and the one-dimensional Schrodinger operator with a Dirac delta 
potential. Furthermore, the same result holds when we impose other boundary 
conditions such as Neumann, periodic, and Floquet boundary condition on the 
longer sides of the rectangle. In Section we review a one-dimensional model 
which turns out to be closely related to low-energy eigenfunctions of Seba billiards. 
In Sectionj^ we consider the optimal estimate measuring the ratio of the localization 
induced by a point scatterer with an arbitrary parameter. Then in Section we 
present the main theorem by combining the results of two previous sections. In 
Section we will give numerical results regarding the asymptotic behavior of the 
localization mentioned in the main theorem for the Seba billiard. Rigorous proofs 
of the main theorem and the intermediate steps are attached in Section]^ 


2. One-dimensional model with a Dirac delta potential 


Let 0 < xo < a, c S M and consider a Schrodinger operator —A — cSxo on [0,a] 
with the Dirichlet boundary condition. We assume that ^ is irrational so that 
(t(—A), the spectrum of the Dirichlet Laplacian, and cr(—A — cS^g) are disjoint. 
If ^ is rational, then there exist Laplacian eigenfunctions sin (^) vanishing at 
Xq for some integer n > 1. Since these eigenfunctions do not feel the presence of 
the delta impurity, they remain as the eigenfunctions of —A — cSxg with the same 
eigenvalues. If ^ is irrational, then sin does not vanish at Xq for any integer 
n > 1 so it cannot be the eigenfunction of —A — cS^g ■ The eigenvalue has to change 
accordingly as well. This implies that 

cr(—A) n (t(— A — cSxg) 0 if and only if — G Q. 

In other words, ^ should be irrational in order to keep cr(—A) ntT(—A — cSxg) = 0. 
The eigenvalues z G a{—A — cSxg) are given as the solutions to 

(2.1) yz (cot (v^xo)-I-cot (^/^(a — Xo))) = c 


having multiplicity 1 and the corresponding eigenfunctions G L^([0,a]) 

( 2 . 2 ) 




N [cot(-\/zxo) sin(-y/z(x — xq)) -I- cos(-\/z(x — Xq))] , 0 < x < xq 

N [— cot{^/z{a — Xo)) sin(-y/z(x — xo)) -I- cos{^/z{x — Xq))] , xq < x < a 


where A > 0 is the normalization constant so that Hiaqo.a]) = 1- Note that the 
superscript (1) of indicates that the eigenfunction is from the one-dimensional 
model. For each c G K, we denote the eigenvalues of —A — cSxg by zi^c < -^ 2,0 < 


■^3,c *\ * * * ■ 

Now we consider the limits of eigenvalues and eigenfunction as c 


00 . By (2.1), 


we obtain an increasing sequence z„ 
two distinct sets, namely. 


= lim, 


c—>-oo '^n,C5 _ 


n > 1 defined as a union of 


(2.3) 


{-^ 1,00 5 ^ 2,00 J ^ 3,00 5 ' ' ' } *^1 ^ *^2 
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Figure 1. A schematic graph of LHS of (2.1) as a function of z. 
Zn,co in (2.3) are marked as dashed vertical lines. 


where 


(2.4) 




m > 1 


5,= 


a — xq 


m > 1 


Remark. For a fixed c G K, (13 provides an increasing sequence of eigenvalues 
Zn,c >0; n = 1, 2, • • • .. On the other hand, for any z G (0, oo) \ {z:„,oo I n > 1}, 
there exists a unique c € M also given by (2.1) such that z G cr(—A — cS^g). See 
Fig. a In addition, since 0 < ^ < 1 is irrational, {zn,oo | ^ > 1} interlaces with a 


sequence i l^j 


{(^)^|n>0},i.e., 


(2.5) 


0 OO ^ -^2,00 ^ -^3,00 ^ 


As c —>■ oo, the eigenfunctions tend to localize either on [0, a;o] or on [xq, a] in the 
L^-sense. We now allow c to take the value +oo and let V'ilh (^) = lin3c->-oo {^)- 
More precisely, if z„_oo G Si, then 

(x) = lim 4^') (x) = I yS - ^o)) ’ 0 < X < xo 

c^oo I Xo < X < a 


If Zn,oo G S 2 , then 


(x) = lim (x) = 

n,<x> C—>-00 

Hence, 

( 2 . 6 ) 


0, 0 < X < Xo 

, {■\/^n,oo{x - Xo)) , Xo < X < a 


lldi!oolU^([o,a^o]) _ 11. if ^n,oo ^ 


IldlhlliAIO.a]) 

On the other hand, for z > 0, has the L^-expansion: 

sin(^xo) 


0? if ^n.oo ^ ^2 


(2.7) 


= M^CnSm(J^x'^ , Cn = 


(w) - 


n— 1 " \ ^ I Z 

where M is the L^-normalization constant. Note that if we assume that — and 

a 

za? are constant, then 

(2.8) M oc a-^/2 
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Since 


l=ll^i^^llL[0.a]=^^E h 


^ / a:o V 

sin ^nTT— J 


n—l 


UTT 


C 


“ = CM^a^ 
2 


for some constant C. We will use (|2.8[) to prove Lemmain Section]^ 


3. A POINT SCATTERER ON A RECTANGLE 


3.1. Boundary conditions. We consider a point scatterer with various boundary 
conditions for which the presence of a point scatterer induces the same type of 
localization of eigenfunctions as in [ 8 ] so that the Seba billiard becomes one example 
of this generalization. 

More precisely, we continue imposing the Dirichlet boundary conditions on [0, a] x 
{0, b} C dn and assume an arbitrary boundary condition on {0, o} x [0, b] C d^l 
satisfying the following properties: For — A on with the boundary condition, 

the eigenvalues Ai < A 2 < A 3 < • • • and the eigenfunctions - ‘ ‘ G 

are equal to 

(3.1) Xn = Xu„u, = n>l,ni>l,n 2 >l 

and 


(3.2) (l>n{x,y) = (I>nun 2 {x,y) = sin (^ 2 ;) gn^ , 

0 < X < a, 0 < y < b, n > 1, ni > 1, n 2 > 1 

for some 0 < < P 2 A • • • and a set of nonzero orthogonal functions | n > 

1} C L^([0,1]) such that 

(3.3) #{Pn \ n > 1} = 0{y/v), V ^ OO. 

For example, the Dirichlet, Neumann, periodic and the Floquet boundary con¬ 
ditions on [0,a] X {0,5} C cAl all satisfy these conditions. More precisely, 

Au + Xu = 0 in n, A € K, u{0, y) = u{a, y) = 0, 0 < y < 5 

has the eigenvalues A„ and the eigenfunctions (/>„ determined by and y„ as in 

(1) Dirichlet: u{x, 0) = u{x, 5) = 0, 0 < a: < a 

i/„ = {jnr f , y„(y) = sin(n 7 ry), n = 1, 2, 3, • • • 

(2) Neumann: dyu{x, 0) = dyu{x, 5) = 0, 0 < a: < a 


(3.1), (|3.2[) for each boundary condition on [0, a] x {0,5} C dXl. 




{mrf, gn{y) 


cos(n7ry). 


n = 0 

n= 1,2 ,--- 


(3) Periodic: it(x,0) = u{x,b), 0 < x < a 

Vn = i2mrf , gn{y) = n e Z 

(4) Floquet: u{x,0) = e“*®u(x,5), 0 < x < a for some 9 G (—7r,7r) 

= (2n^ + ef , y„(y) = e*(2-+«)«, „ g Z 

For (2),(3) and (4), we define I'm 9n (n > 1) by rearranging i>„ (n > 0 or n G Z) in 
the nondecreasing order with a new index n > 1 so that (3.3) holds. 
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Figure 2. Geometry of a point scatterer at Xq = {xo,yo) (marked 
as x) in 17. The left part of the plate divided by Xn is denoted by 
17i = [0 ,xo]x[0,6]. 


3.2. Spectral properties of a point scatterer. Now we construct a point scat¬ 
terer at Xq = {xo, |) on a rectangular plate 17 = [0, a] x [0, &] as in Fig. with 
the boundary conditions discussed in Section [3^ We continue assuming ^ to be 
irrational as in Section [2] 

First, consider the Laplacian — A on L^(17) with the boundary condition satis¬ 


fying (3.1), (3.2), (3.3). Then restrict its domain to the functions vanishing at Xq. 


By the theory of self-adjoint extension developed by Von Neumann, such a sym¬ 
metric operator has a family of self-adjoint extensions —Axo,n,a with a parameter 
a € (— 00 , 00 ]. The parameter a is following Albeverio’s notation from |10j as in 
our previous paper [8]. Note that a determines the strength of the point scatterer 
at Xq € 17 although it is not equal to the coefficient c € M of the delta potential 
in More precisely, — A -|- cJxq needs renormalization to be properly defined 

as a self-adjoint operator on L^(17). Let Gz be the integral kernel of the resolvent 
(—A — z)~^ : L^(yt) —>■ namely. 


Gz{x,x') = 


^n(x)(/)„(x') 

t ^ A — z 

n=\ ^ 


so that for / S T^(17), 


(-A-z) V(x) = [ G^(x,x')/(x')dx'. 
Jn 


Then the integral kernel of — A -|- cJxo; c G K formally reads 


(-A-I-c<5xo - z) (x,x') = G 2 (x,x')- 


- - G^(xo,xo) 


1 -1 


G2(xo,x')G^(x,xo). 


However, this approach fails if G 2 (xo,Xo) diverges which did not occur in the one¬ 
dimensional case. Therefore, we have to discard c and renormalize the delta po¬ 
tential with a new parameter a € (— 00 , 00 ] to obtain the resolvent formula for the 
point scatterer —Axo,n,Q. For z G C \ cr(—Axo,n,a), the integral kernel of 

(~Axo,n,a ~ ^) ^ • T^(17) —>■ L^(fi) 

is given by 

(3.4) (-Axo,n.a - 2)"^(x,x') 
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Figure 3. A schematic graph of F{z) of (3.5). A,, 
marked as dotted vertical lines. 


in (3.1) are 


= G^(x,x') + [a - F(z)] ^G 2 (xo,x')G 2 (x,Xo), ae(-oo,oo] 


where 


(3.5) 


Fiz) = 

n—1 


1 


A. 


Xn — Z A^ + 1 


is derived by applying Von Neumann’s self-adjoint extension theory [5] to the Lapla- 
cian restricted to the smooth functions vanishing at Xq. 

According to this parametrization, the point scatterer disappears in the norm 
resolvent sense of becoming — A as a —)■ ±oo whereas it gets stronger when laj <C 
00 . In particular, a = oo corresponds to the Laplacian with the same boundary 
condition on 517 so we will only consider a G K for which the presence of the point 
scatterer actually perturbs the system. We may also interpret a as the inverse of 
the coupling constant vb or vg in Shigehara’s work HE]- For more details on the 
construction of (3.4), (3.5), see m- 

Now consider the eigenvalues of the point scatterer with the boundary conditions 
in Section 3.1 For simplicity let us omit Xq and 17 from —AxQ,n,ct, namely, 


-A„ = -A 


Xq ,n,Q; 


since those quantities are fixed. Also, we denote mult(z, P) by the multiplicity of 
z as an eigenvalue of an operator P. For a G K, we can divide ct(—Aq.) into two 
types: 

(1) Perturbed eigenvalues: cr(—A q.) \ cr(—A) and 

(2) Unperturbed eigenvalues: a{—Aa) n cr(—A) 

where each of them is obtained by different conditions. For a G K, 


(3.6) z G ct(—A c) \ cr(—A) if and only if a = F{z) 

with F{z) defined in (3.5). See Fig.Furthermore, mult(z, — Aq.) = 1. 
On the other hand, define fj, and /Iq 


fj,(z) = mult(z, —A) = #{n > 1 | z = A„} 

= #{n > 1 I z: = A„,())„(xo) = 0}. 

Then for a G K, 

(3.7) z G cr(—A q) n cr(—A) if and only if ^J^o{z) > 1 or ^(z) > 2. 
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with the corresponding eigenspace 


(3.8) 


^ ^ Cn4^n 


^ Z — \n 


Cn0n(Xo) =0, Cn GC 


Z — Xn 


By (3.71 and (3.8 1 , we observe that the unperturbed eigenvalues and the corre¬ 


sponding eigenfunctions both are independent of a S K. The proofs of (3.71 and 


(3.8) can be found in Chapter 2 of m with generalized statements for a point 
scatterer on a compact Riemannian manifold of dimension two or three. 

Let < Z 2 Q, < < • • • be the eigenvalues of —Aq,. Now we remark some 

properties of cr(—A q). 

Remark. For any z € IR\ {A„ | n > 1}, there exists a unique a G K given by (3.6) 
such that z G tr)—A q). For any a G K, cr(—A q) interlaces with (t(—A), namely, 


(3.9) 




< Ai < z^_Q < A 2 < 


In addition, for n > 1, we have 


lim z(j Q = Xn, 


lim 


/ 

^n-\-l,oc 


— Ar] 


lim 


4,0. = 


/Q\ 

Let tpz be the eigenfunction associated with z G ct(—AQ)\ cr(—A). Note that the 
superscript ( 2 ) of ipz ' indicates that the eigenfunction is for the two-dimensional 
model. Then we have 

= MGzixo,yo; x,y) 

where is the integral kernel of the resolvent operator (—A — z)~^ : L^(n) —>■ 
L^(n) and M > 0 is the L^-normalization constant. In addition, has the 
L^-expansion up to a multiplicative constant: 


(3.10) 


= E 


n—1 


yo) 

Arj, Z 


= E E 

m—l n2—l 


4^711,n2 (^Oj yo) 

^rj.i .r).o ^ 


'■^TLi ,712 


Note that for a G K, any / in the domain of —Aq has a logarithmic singularity at 
ixo,yo) or vanish at (a;o,yo)- More precisely, as {x,y) -)■ {xo,yo), 

(3.11) f{x, y) = Cl (^In [x - Xq)"^ + {y - yoY + ^ 2 ( 0 ?)) + o(l) 

for some Gi G K and C2{a) G K fixed for each a. Proofs can be found in Chapter 1 

of m with a general statement for a complete and smooth Riemannian manifold 

( 2 ) 

of dimension two. This implies that ipz also contains the logarithmic singularity 
unless ipf'\x,y) = o(l) as {x,y) (a;o,?/o)- 

Remark. Without loss of generality, we may assume area(r2) = a6 = 1 or equiv¬ 
alently, 

a = y/E, b = 

Ve 

(‘2'\ 

by scaling fl, z(j q and ipp/ simultaneously according to the following properties: 

Consider point scatterers — Axo,n,Q and —Arxo,rQ,a on 11 = [0,a] x [0, &] and ril = 
[0,ra] X [0,r6], respectively, with a coupling constant a. Then 

(3.12) z G cr(-Arxo,rn,a) if and only if r^z G cr(-Axo,n,Q_/3) 
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where /3 is a constant defined as 




"”(2) (2) 

In addition, let ipz ’ and be the eigenfunctions corresponding to the eigenvalues 
z £ — and r^z £ — Aq Q,_,g, respectively. Then for some C ^ 0, we have 


= C-ip 


( 2 ) 



Now we estimate the localization of eigenfunctions over fli = [0, xq] x [0, 6] in the 
L^-sense where the perturbed eigenvalue is bounded below by the lowest eigenvalue 
of the Laplacian and bounded above by min{Ai „ | Ai,„ > Ai i, n > 1}. See 
Section for the proofs. 


Lemma 1. Define A 


ni ,712 ’ n2 


as in (3.11 and let a = '/E, b = 


Let 


(3.13) A = min{Ai,„ | Ai,„ > Ai,i, n > 1}, n = min{n > 1 | Ai^„ > Ai,i}. 


Let z £ (Ai,i,A) n cr(—A(j) \ (t(—A) and assume z = z(E) = ^ + v^E for some 
constant c £ K. Then 




( 1 ) 


(3.14) 

and 

(3.15) 


O 


ll'^z ^llL^(n\ai) _ IIV'2-t/,^llL2([a:o,a]) 


( 1 ) 


o 


E^^+fi^CoE-z^ 


^E^^+fi^CoE-z^ 


as E ^ oo, where 
(3.16) 


Co = sup ^, 

n>l 


with £ L^([0, v^]) of (2.21. Note that z and its bounds Aiu,A all depend on 

E. 


4. Main Theorem 


In this section, we provide a sequence of a’s and eigenvalues inducing the local¬ 
ization of certain eigenfunctions by combining (2.6) with Lemma 


Theorem 2. Let 


Ne = 


E — El 


c;2 -p 1 


where v = min{i^„ | > i/i, n > 1}. Choose E sufficiently large so that Ne > fi- 

Define Si, S 2 and Zn,oo as in { 2.4) with a = '/E. For n = h, - ■ ■ , Ne, there exists 
a unique parameter a„ explicitly defined as 


an = ^ \(t>n’{xo,yo)\'^" 


Xn 


Xn' l^\E Zji—1.00 


A£, + l 


(4.1) 
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such that 


^n,an — 1,00 ^ ^a„)- 


In addition, as E ^ oo, 
(4.2) 

llL2(n\ai) 


ll<^„ IU=(0) 


= o 


1 


= o 


y/E^{fi^Co - Ui) - i;7r2(n2 _ i)^ 

'_1_) 

^/E^ifi^Co - lyi) - E7T^{n‘^ - 1 ) J ’ 


^n^oo ^ * 5*1 




Note that h in ( |3.13 ) is equal to 2 if z/i ^ which is the generic case. Then 
Ne > h is chosen as the largest number satisfying 


, , /mr\'^ vi 


and 


<t>n{x,y) = sin 


\ a 


’ \ 

-j 5i(y): 0<a;<a, 0<y <b 


for all n < Ne- Since ^ and gi are fixed with respect to n, we can say that Ne 
describes how many low-level modes of the 2D harmonic vibration on D can be 
approximated by those of the ID harmonic vibration on [0, a] with some constant 
quantities such as ^ and gi. By Theorem this idea can be extended to approxi¬ 
mating a point scatterer on D by a Schrodmger operator with a delta potential on 
[0,a] in which the localization of eigenfunction has been already exhibited. How¬ 
ever, the lowest mode of the point scatterer does not have a corresponding mode 
in the ID model so n < h (or n = 1 in generic cases) should be excluded from 
the argument. See the proof in Section for more details. Note that localized 
eigenfunctions could appear sporadically above Ne as presented in [5]. However, 
such cases will eventually disappear as n —>■ oo since the high-energy eigenfunctions 
of Seba billiards tend to localize rather in the momentum space. See [T^ for the 
details. 


Remark. In the asymptotic notation (4.2), n may be a variable depending on E. 
If n > 2 and n = o{E) as E ^ oo, then (4.2) reads as E ^ oo. 


(4.3) 


l|L2(o\ni) 




= O 


{E-i) 


if ^n,oo ^ ‘5*1 


= o 


{E-i). 


if 


€ 52. 


11^2(0) 


In addition, since 
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n 


Figure 4. Plot of 2 < n < Ne given by (4.4) where Xg = 
E = IOtt with the Dirichlet boundary condition on dfi. 


(4.2) actually implies that there exists a constant C > 0 and M > 0 such that for 
all E> M, 


n < ^ ^ f 1 

■ ^ [vE^{fi^Co-i^i)-E7T^n^-l) 

||V'if^^llL2(n) \^yE^{n‘^Co - ui) - Eir'^^n^ - 1) J 


if ^n,oo ^ *5*1 


if ^n,OD ^ *5*2 • 


In particular, we may apply Theoremj^to the examples mentioned in Section [3.1| 
and calculate a„, z!^ that generate the localization of modes as exhibited partly 
in [S]. Note that the Dirichlet boundary condition is given on {0,a} x [0,6] for all 
cases. 


Example 1 (Seba billiards). For the Dirichlet boundary condition on dVt, 

(4,4) a„= f; 5: ^ 

m—l n2>l 
712 odd 

and 


*^ni,n2 7t‘^E Zn^cfO ^ni,n2 ^ 


where 


+-^n-l.oo G 0 '(-Aa„), 2<n< \/ 3E'^ + 1 


/ f* 1 /( Q O 

Ani.na = + ni-K E, Til >1, 712 > 1- 


See Fig. 1^ for the plot of a„ given by (4.4) . 

Example 2. For the Neumann boundary condition on [0,a] x {0,6} C i9D, 



712 even 


= Z„_yoo G cr(-Aa„), 2<n< \/i?2 + 


and 
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where 




A. 


— E ^1 — ^5 ^2 ^ 0 - 


Example 3. For the periodic boundary condition on [0, o] x {0, b} C dfl, 

(4.6) a„ = y y sin^ (^- , ] 

ni = ln2GZ \V^ / \A„i,„ 2 - Z„,oo ^ni.ns+l/ 

and 

4.a„ = 2 :n-l,oo e ct(-A„^), 2<n< y/4£’2 + 1 

where 

2 2 

Ani.n2 = + 471271^F;, Tli >1, 712 € Z. 

h/ 

Example 4. For the Floquet boundary condition on [0, o] x {0, b} C d^l with some 
9 £ (-7r,7r), 

^n' \ 



m—i 712 


'n,oo H“ 1 


and 


where 


4,a„ = + Z„_yoo G CT(-Aa„), 2<n< 


/4£;2 




+ 1 


Am,712 “ y=;-^ (271712 + Til >1, 712 G ^• 

E 

5. Numerical Results for Seba billiards 

In [8], we numerically exhibited several kinds of localized eigenfunctions of the 
Seba billiard without specifying a„, the parameters localizing the n-th mode up to 
an error of 0{E~^) to one side of the rectangle. We now numerically check how 
accurately the localization occurs for —Aa„ for E large as expected by Theorem]^ 
For the sake of convenience, let the x-coordinate of the point scatterer be fixed at 
xo = y. Note that, however, the qualitative properties we observe also hold for 
other values of xq. 

First, consider a Seba billiard with a fixed eccentricity E = IOtt. Although a S M 
is the variable to be considered for —Aq,, we can let the eigenvalue z G cr(—A^) 
itself be an independent variable inIR\{A„ | n > 1} and let a G K depend 
on z since for each z, there exists a unique parameter a G M defined by (3.6) 
such that z G tT(—A^). In this point of view. Fig. ^ shows the L^-norm of the 
eigenfunction i/jz ' over fli as a function of z. Note that we assumed \\tjjz \\L^(n) = 1 

so liy^ 11^2(01) = I and ||y^||L 2 (Q^) = 0 mean that gets completely localized 
in fli and fl \ Oi, respectively. The dashed and dotted lines indicate the sets 
of optimal eigenvalues Si + tt'^E and S 2 + exhibiting the localization in fli 
and n \ Oi expected by Theorem One can observe that those lines in Fig. 
thoroughly estimate the n-th lowest eigenvalues (n = 2, - • • ,Ne) of the localized 
eigenfunctions. For example. Fig. shows the partial sum of the first 10® terms of 

(3.10) where z = Zg^^^ = tt'^E + Z 2 ,oo. Since Z 2 ,oo = ^ eigenfunction 
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Figure 5. (Seba billiard) as a function of z where 

xo = ^, E = IOtt. The dashed and dotted vertical lines indicate 
sets of eigenvalues: Si + tt'^E and S 2 + •k'^E, respectively, in which 
the localization expected by Theorem [2] occurs. 



Figure 6. (Seba billiard) Plot of computed as the partial 
sum of the first 10® terms of (3.10) where E = IOtt, Xq = ^ 
10^. 


1.78, z = z'ii 


3,123 


3.13 ■ 


tends to localize on fli. A small peak at (ccq, |) indicates the logarithmic divergence 
of mentioned in ( |3.11 ). This phenomenon was not emphasized in the figures 
of localized modes in [S] . 

Second, we introduce numerical examples verifying how strongly the eigenfunc- 
tions Ip]., localize in terms of the asymptotic estimate in (4.2) as if —>■ 00 . For 
E > 0 and for n = 2, • • • , Ne, consider an operator — Axo,n,a„ and its n-th lowest 
eigenvalue where Xq = (xq, |), 11 = [ 0 ,if 2 ] x [0, if”^] and is given by 

(|4.1|). We choose xq = 0.3-, 0.7-, 1.1-, 1.5- to observe how the location of a 
point scatterer affects localization of eigenfunctions. (See Fig. 0) Define Cn.E as 


l|L2(n\ni) 


lU^(ni) 






.,( 2 ) 


IU2(0) 


if -2^71,.00 ^ 


if ^ 71.-00 ^ So 















































LOCALIZATION OF SEBA BILLIARD EIGENFUNCTIONS 


13 




(a) xo =0.3f « 0.10-a 


(b) xo = 0.7f « 0.22 -a 




(c) Xo = l.lf ~ 0.35- a 


(d) Xo = 1.5f « 0.48 • a 


Figure 7. (Seba billiard) Rate of localization of eigenfunctions 
; (2 < n < y/SE'^ + 1) approximated with data near E = 
2tt, 4?!, Gtt, Stt, IOtt with various xq. Eigenfunctions localizing in 
and n\fli were marked as x and •, respectively. The dashed lines 
correspond to the theoretical bound k = —1.5 as E —> oo given by 
Theorem 


and define the rate of localization as a number k such that Cn.E = 0{E^) as E —>■ oo. 
Note that /c < 0 for some n>2 implies that the L^-norm fraction of the localized 


61 decays 


eigenfunction '0^' q unlocalized region (e.g., the right part of Fig. 

^ ^ ( 2 ) 
as E increases. In addition, |fc| measures how fast the unlocalized fraction of ip),/^ 

diminishes as E —>■ oo as long as fc < 0. 

The numerical results indicate that fc < — | for any n > 2 as predicted by 
Theorem I Fig. shows the rate of localization of as a function of n with 

various xq. Data points for eigenfunctions localizing in fli and D \ Oi were marked 
as X and •, respectively. According to (4.31, all points in Fig. [^should lie on 
or below the horizontal line of — | no matter where the localization occurs in the 
rectangle. In particular, one can numerically observe that the rate of localization 
converges to the theoretical bound marked as the dashed line as E —>■ oo 

Note that keeping the irrational ratio between xq and a as assumed in Section 
is crucial to localize all eigenfunctions as E ^ oo. For example, placing a point 
scatterer at the midpoint ccq = | does not induce localization of eigenfunctions 
since fljand D \ are symmetric to each other. However, one can observe that 
in Fig. 7d all eigenfunctions localize within a controlled error better than 0{E~'i) 
although *0 = 1-5 “ ~ 0-48 • a is close to the midpoint of the rectangle. This implies 
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that even a slight change of xq can induce a drastic change of eigenfunctions in 
terms of the localization ratio. 


6. Proofs 


Proof of Lemma\^ We rewrite (3.10) using (2.7) as 


OO OO 


^ ^ (l)m,nAx,y) = ilJiov,{x,y) + •0high(a;,i/) 


+ l^noE - z 


where 


ni = l 712 = 1 E ' ^n2 


V'iow(a;,2/) = X! X! 


^ni ,712 (^0i Vo) 


m—l7l2 = l E ^ 


’^71i,7l2 {^1 y); 


OO OO 


/ / ^ 0711,712(^0, yo) , / X 

0high(3^,y)— / ^ ^ 7r2r,2 „ 0711,712 {^,y)- 


ni = l7i2=n £; “r ^712 


+ Un.E - Z 


By (3.21, 0ni,7i2 can be rewritten as 


0711,712 (x, y) = sin 


In addition, by (2.7), 


OO n— 1 


fjlo^{x,y)= Ki,U 2 {x,y) 


T^Tl^ 

711 = 1 7l2 = l ^ ^712 

n —1 


+ r-njif - 2 
1 


9n2 (yEyo)jjf^i-^^^Ei^)9n2 iyEy) 

2 = 1 

(ill 5^2 (v^yo)gn2 (v^y) 


V 712 = 1 


Note that VEyo = ^ is a fixed quantity for all if > 0 and gi,g 2 ' ■ ■ given in 
( |^ are orthogonal. By applying ( [^ to obtain ||V' 1 ow||l 2 (o) = CE 

for some constant C > 0 since \\yi^^E\\L'^([o,a\) = 1 and 


n—1 


1 I ^ 

IIV'low||i2(0) = J^Wy^EWhiXOM) XI |y"2(v^yo) 5 n 2 (v^ •) 


( 6 . 1 ) 


712 = 1 
2 


L2([o,6]) 


= X |y"2(v^yi 

712 = 1 

= C^if^ for some C > 0. 


yn2llL2([0,l])-' 


\E ^ ] for some c 


On the other hand, by (3.3) and z < Ai^ 2 , there exists C > 0 such that 

2 


II / l|2 I \(l>ni,n2{xO,yo)\ \ n , 

IIV’high|L2(n) = 2 ^ 2 ^ -^2 --- ll'/'ni, 


7ii = l 712—n \ E ' ^n2^ ^ / 


n2llL2(a) 


OO OO 


^^X X 


711 = 1 7i2=n \ E \ >^2^^ 


10711,712(^05 yo)| 

^+rilCoE-z, 


) 2 

nOO nOO 

<C / + 

J^JnVGdE 
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< c 


'o J^J^+n^CoE (r2 - z) 


-—= C , ^ —- 

2 ^ + h'^CoE - z 


where h, Co are given by (3.13) and (3.16). Therefore, 

1 


( 6 . 2 ) 


||'*/'high||L2(n) < C 


^ + h^CoE - z 


By ( |6.1[ ), ( |6.2[ ), we obtain 

IIV’z ^ 11^2(121) ^ ll'*/’low||L2(ni) + ||V'high||L2(ni) 


11^2(0) ll'01ow||L2(f2) IIV^highlli^JQ) 

'IIV^high||L2(ni) 


||i/^1ow||l2^Q^ 

^ ||'01ow||L2(f2^) ^ ^ 

~ ||'01ow||l2(q) 

^ \\'>I^''z-uie\\lH[0,xo]) 

Hi^-^,EhH[0,a]) 


||'^iow||L2(n) 

1 



E^I%+n^CoE-z, 


Ej^+fi^CoE-z, 


and similarly, 


ll'0z^^l|L2(nj) ll^i-^iulU^dO.xo]) 




^ II^W 


-c 


lL2(n) 


^_,.,i5llL2([o,a]) 


Ej^+n-^CoE-z, 


for some C > 0 as C —)■ cxd. This concludes the proof of (3.14). Similarly, we can 
prove (3.15) by switching fli and ft \ fli. □ 

Proof of Theorem^ Note that a = '/E and b = since area(n) = 1. Then we 
have A„ = A„,i for all n < Ne- In addition, (3.9) implies that for n = n, • • • ,Ne 
and for any a S K, satisfies 

Ait = Ai < < Aat^, = \ne,i < '^1,2 

so we can apply Lemmaj^to z'^ ,^ S tT(—A q.) \ (t(—A). Furthermore, for each n, we 
may choose a specific a = a„ so that the localization of is maximized. More 

preciesly, consider the 1-dimensional model discussed in Section with a = \/~E. 
By (2.5), we have 

Ait < ^1,00 + viE < A2T < ^2,00 + viE < Ast < ^3,00 + viE < • • • 


Hence, for each n = n, - ■ ■ ,Ne, there exists a unique a„ S M given by (4.1) such 
that 


^n,a„ — Zn-1,00 + EiE e (Si U S2) + ViE 





























and by (2.6), 


' II 

-vieWl'^UxoM) 

= 0 

if 4, 

ll^<„. 

^-^^E\\L^[0,a]) 

II 

-vie\\l'^{[0,xo\) 

= 0 

if 4, 

. ii'^Sc 

-mE L2([0.a]) 

1 ^n,a 

^ \ri -I ^ ^ 

n < Ne and 

An-l = 

7r^(n — I)^ 

Xn — 


E 


viE. 


Therefore, we can apply Lemmaas follows: If — viE G S 2 , as E ^ oo, 


WL^iili) IIV'z'^ -ly^EWL^ilO^xo]) 




( 2 ) 


lL=(n) 




o 


vxE\\L'^{[0,a\) 


E 


= 0 + 0 


y y^E^{n'^Co — vi) — E'K‘^{n’^ — 1) ^ 


If. 


— eiE g Si, as E ^ 00 , 


ll^2^^JU+n) ||V'i;^^^_^i£;llL+[ 0 ,a]) +n^CoE- z: 


+ 0 


/ 

n,a.ri 


= 0 + 0 


y/E3{fi^Co - El) - En^in^ - 1 ) 


□ 
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